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Nanosized spherical system of magnetic moments interacting indirectly via the RKKY mech- 
anism is studied. The interaction energy that determines the temperature of the ferromagnetic 
ordering, depends strongly on the system size. Obtained in the mean-field approximation, dimen- 
sional and concentration dependencies of the Curie temperature testify to the necessity of taking 
into account the finite size of such systems to calculate their features. Results may concern both 
artificially constructed nanosystems and naturally arising formations (such as clusters of magnetic 
ions in diluted magnetic semiconductors, etc.). 



In systems with free carriers of high concentration (metals or degenerate semiconductors), 
indirect magnetic impurities' interaction of Ruderman-Kittel-Kasuya-Yosida (RKKY) type is 
considered as one of the basic mechanisms of the magnetic ordering [T]. There are a lot of 
papers dealing with RKKY interaction in three-, two- and one-dimensional systems of infinite 
size in El El SI • However, to our knowledge, nobody considered how that interaction should be 
modified for systems of the finite size. In the present paper, we consider this problem for the 
case of the spherical system of the finite radius. 

Estimating the energy w{r-) = — J(r)SiS2 of the indirect RKKY-interaction of magnetic 
impurities with spins Si, S2 spaced at the distance r is based on making use of the expression 



k-l-q 

obtained in the second order of the perturbation theory [T]. Here, is the number of lattice 
sites, Jpd is the exchange energy for the interaction of the impurity spin with a free charge 
carrier, is the carrier energy, /(-Ek) is the Fermi-Dirac function which in the degenerate 
case equals /(-Ek) = I at k < kp and f{Ek) = at k < kp, where kp is the Fermi momentum. 
The prime by the the sum over q means that q 0. The exponent e"'"/^ in ^ reflects the finite 
carrier mean free path I. 

In the continual approximation, the summation in (|T} is replaced by the integration which 
performed usually over all k < kp and |k -|- q| > kp. Then, the standard result corresponding 
to the case of the infinite system reads 

J(r) = — /o$(r) exp(— r//), $(r) = ^-j [ip{r) cos ip{r) — simp (r)], (2) 

where 

^ = ^(^^i)' Vir) = 2kp,r, (3) 

a is the lattice constant, kpQ = (Svr^p)^/^ is the Fermi momentum of carriers of the concentra- 
tion p. 
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In the case of finite system sizes or magnetic ions' clustering, the classic expression Q for the 
energy of the RKKY-interaction should be rectified. For simplicity, we consider the case when 
magnetic ions form the spherical cluster. Due to the quasi-neutrality, its radius R determines 
not only the area where ions are arranged but also the region where carriers, produced by 
those ions, are localized. In other words, the carriers are contained in the potential well of the 
radius R. Therefore, the carrier momentum k and its variation q are limited by the intervals 

h < k < kp, ki < q < k2, (4) 

where 

ki ^ vr/i?, k2 ~ vr/a. (5) 

In addition, due to the spatial quantization the distance between energy levels of carriers grows 
that leads to increasing the Fermi energy and Fermi momentum with decreasing the well size: 
kp = kp^R)- Together, it complicates calculations and the final expression for J(r) turns out 
to be more bulky than the canonical expression 

The finite mean free path / results in smearing energy levels of carriers due to their collisions. 
Therefore, the lowest value k is defined by the system size and equals ki ^ n/R only if the 
collision broadening h/r of levels is less then the energy n'^h^ /2mR^ of the first level. That 
condition could be written in the form 



2k 



i?2 / 



^ > (6) 



meaning that our approach relates to the small enough systems only. If, for instance, //a = 10, 
akp = 1 then R < 10a. In the general case, one could use the value 



as the left boundary of inequalities ^ . 

To proceed one should estimate how the Fermi momentum depends on the cluster size. 
The total number of free carriers in the cube of the size R with the spherical Fermi surface 
is defined by the number of cells of the volume (27r)^ in the phase space and in the limit of 
kpR — > oo equals Np ~ k\R^ / {2'kY . For a finite A;i?i?-value, the number of carriers is defined 
by the number of points in the wavenumber space with coordinates divisible by {2n/R). As 
none of those coordinates could equal zero, in that case 

Np ~ {kpR/27rf - a{kpR/27rf, 

where the correction (proportional to a ~ 1) is associated with the "forbidden" points, posi- 
tioned in the coordinate planes of the wavenumber space and with the "correct" points neigh- 
boring to the Fermi surface but not fallen inwards, as well. For the carrier concentration 
n = Np/R^, it follows n ~ kp{l — a/kpR). For n = const and kpR ^ 1 that leads to 

a 

kp{R) - kpo ^ — , 

where kpQ = kp{oo) is the Fermi momentum in the infinitely large system. Everywhere below, 
we use kp to mean the value 

kpiR) = kpo + ^. (8) 
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Let us turn now to calculating the energy of the RKKY interaction. Assuming = h'^k^/2m, 
£'k+q = h?{k + qY /2m and designating the angle between r and q as a, and the angle between 
k and q as 6*, we obtain 



J(r) 



exp(— r/Z) — 



iqr cos a 



E 



f{Et 



2k cos 9 + q 



or, in the continual approximation 



k2 n kp IT 

J(r) = -/oexp(-r//)a^/" qdq I e"^''"''' sin a da [ k^dk [ . 
TT J J J J 2kcos9 + q 

fci fei 



(9) 



(10) 



Non-complicated but laborious calculations lead to rather cumbersome result which could 
be represented in the relatively simple form by the help of the operator L determining the value 
of the double definite integral J^^ dk J^^{...)dq with the primitive function g): 

L^{k, q) = ^(A;i, A;i) - -^{kp, k,) + ^{kp, k^) -^{k^, k^). 



Then 



where 



J(r) = - Joexp(-r//) {akpYlm r{k,q), 

71 
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ci{k, q) - ^2{k) C2{k, q) + $3(fc, q)] 



{2kFrY 

$i(fc) = 2kr cos2kr — sin2/cr, ^2{k) = 2kr sin 2kr + cos 2kr, 



(11) 

(12) 
(13) 



ci(fc, q) = Si[{2k + q)r] - Si[(2A; - q)r], C2{k, q) = Ci[{2k + q)r] - Ci{\2k - q\r), (14) 



$3(fc, q) = cos qr{l + 2A; r ) In 



2k + q 



2k -q 



+ gr (sin qr — —qr cos qr) [ In 



2k + q 



2k -q 



Ak 



(15) 



Si(x) 



sint 



dt, Ciix) = — 



cost 



dt. 



(16) 



To the traditional situation {R — > oo) there correspond ki — > 0, A;2 ^ oo. In that case Ci vr, 
C2 ^ 0, $3 — > 0, and 



TT ^ ^ 2kFr cos2kFr — sm2kFr 
^i{kF) = TC 



{2kFr) 



{2kFrY 



Hence, (fTT| reduces to the standard expression Q. For finite values fci, ^2, the interaction 
energy J(r) should be calculated with Eqs. (fTI|l -fl5). 

The local effective RKKY- field -^rkky, generated in a given point, is defined by the relation 
yui^RKKY = J2i where rj is the distance from that point to the ^th magnetic impurity, fi is 
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the impurity magnetic moment. In the continual approximation, the sum could be replaced by 
the integral ^H-^kky = f J{r')d?r' where the integration is spread over the volume occupied by 
impurities. Contrary to the infinite system, the value of that integral depends on the position 
of the considered point. For the spherical system, the effective field could be characterized by 
the value 

R 

/^^RKKY = 4vr| J{r)rHr (17) 

of that integral in the center of the sphere. Here is the minimum distance between impu- 
rities determined by the discreteness of the crystal lattice (for instance, the minimum distance 
between extrinsic Mn atoms, replacing Ga atoms in GaAs lattice, amounts r^^^ = a/ \/2 f» 4 A). 

How much would the local field -^rkky be non-uniform inside the sphere one could judge 
by making note that in the case with the interaction length / being comparable or shorter than 
the radius R, the local field at the surface of the sphere should be approximately half as large 
as its value in the center of the sphere. It is not hard to show that for any point inside the 
sphere being offset by the distance h < R from its center, one could employ the relation 

Y r 1, r<R-h 

liH^KKY = ^T^ J{r)r^F{r)dr, F{r) = { _ (^h - rf „ , p^,, (18) 

J y , K h<r < K + ti 

' mm 

instead of (fT7|l . In particular, for the field at the sphere surface {h = R) one could find 

2R^ 

fiHl^^^ = 2nj J{r){l-r/2Rydr. (19) 

''^min 

The results of numerical calculations (see Fig. 1) with r^i^ = a/\/2, / = 3a, kpoa = 1 
show that for R = 10a the part of the sphere where the effective field differs from -f^RXKY 
no more than by 20%, amounts about 85% of its volume, and the average field -^rkky = 
(3/i?3) HKKKY{h)h'^dh ^ 0.93i/^^KKY (the same for R=3a amounts I/rkky ~ l.UH^^j^y). 
Hence, one could, in the first approximation, ignore the non-uniformity of the effective field 
and consider it as nearly uniform and equal to -^rkky ~ -^rkky- 

In the mean-field theory, Curie temperature Tc of the ferromagnetic state arising due to 
the RKKY-interaction is defined by the simple relation 

keTc ~ flHRKKY ~ Ai^RKKY- (20) 

It is of interest to understand how the so-defined Curie temperature depends on the system 
size R at kpo = const or varies with the carrier concentration (determined by the Fermi 
momentum [HI) in systems of fixed (but different) sizes f6]. Corresponding dependencies are 
shown in Fig. 2, 3. 

The dependence Tc{R) turns out to be non-monotone, with pronounced oscillations of the 
period 2a at higher kpo values. They are resulted from that part of the interaction energy J(r) 
which is defined by the function $3(A;,g) given by (15), namely, by the terms singr and cosgr 
with q = k2 = 71 / a. That provides for the observed period. 
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At i? ~ /, the Curie temperature exceeds noticeably its value for i? — > cx) (shown by dot 
lines in Fig. 2). But the most significant distinction of the dependence Tc(-R) (associated with 
accounting the finite ki value) is observed at small R. As soon as the condition © is reached, 
the drop of Tc is observed, as compared to the "standard" (corresponding to ki = 0) dependence 
(shown by dash curves in Fig. 2). That means the impossibility of existing ferromagnetism due 
to RKKY interaction in small clusters. 

As for the dependence Tc{kpo), the Curie temperature increases almost monotonously with 
kpo (i.e., with rising the carrier concentration). In that case also, reaching the condition © 
results in the considerable decreasing Tc (down to zero at some finite size i? ~ /) as compared 
to the Tcikpo) dependence for R ^ oo. 

In conclusion, we studied nanosized spherical systems of magnetic moments interacting 
indirectly via the RKKY mechanism. The interaction energy which determines the temperature 
Tc of ferromagnetic ordering, depends strongly on the system size. Obtained in the mean-field 
approximation, dimensional and concentration dependencies of the Curie temperature testify 
to the necessity of taking into account the finite size of such systems to calculate their features. 
Results may concern both artificially constructed nanosystems and naturally arising formations 
(such as clusters of magnetic ions in diluted magnetic semiconductors [3 El, etc.). 

This work has been supported by Grant No. 06-02-116313 of the Russian Foundation of 
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Captions 



Fig. 1. Spatial distribution of the local effective field -^rkky within the sphere of the radius 
R = 10a and R = 3a at I = 3a, kpoa = 1. 

Fig. 2. Dependencies Tc(-R) of the Curie temperature on the radius of spherical systems with 
various carrier concentrations. Upper panel: I = 3a, lower panel: / = 10a. Dot lines indicate 
the values Tc{R = oo), dash curves show the "standard" (corresponding to ki = 0) behavior of 
the dependencies Tc{R) at small R. 

Fig. 3. Dependencies Tc(fcFo) of the Curie temperature on the carrier Fermi momentum for 
spherical systems of various radii. Upper panel: / = 3a, lower panel: / = 10a. Arrows indicate 
where the condition ^ is reached. 
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